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Abstract 

We analyze the dynamics of neutron beams in interferometry experi- 
ments using quantum dynamical semigroups. We show that these ex- 
periments could provide stringent limits on the non-standard, dissipative 
terms appearing in the extended evolution equations. 



An open quantum system can be modeled in general as being a small subsystem in 
interaction with a suitable large environment. Although the global dynamics of the com- 
pound system is described by unitary transformations generated by the total hamiltonian, 
the effective time evolution of the subsystem usually manifests dissipation and irreversibil- 
ity. 

This reduced dynamics, obtained by eliminating (i.e. by tracing over) the environ- 
ment degrees of freedom, turns out to be in general rather complicated. However, under 
some mild assumptions, that essentially ask for a weak coupling between system and en- 
vironment, one obtains a subdynamics free from memory effects, that can be realized in 
terms of linear maps. Furthermore, this set of transformations possesses very basic and 
fundamental physical properties, like forward in time composition (semigroup property), 
probability conservation, entropy increase and complete positivity. They form a so-called 
quantum dynamical semigroup. [1-3] 

This description is rather universal and can be applied to model a large variety of 
different physical situations, ranging from the study of quantum statistical systems, [1-3] to 
the analysis of dissipative effects in quantum optics, [4-6] to the discussion of the interaction 
of a microsystem with a macroscopic measuring apparatus. [7-9] It has also been proposed 
as an effective description for phenomena leading to loss of quantum coherence, [10-14] 
induced by quantum gravity effects at Planck's scale. [15] The basic idea is that space- 
time should be topologically nontrivial at this scale, manifesting a complicated "foamy" 
structure; as a consequence, transitions from pure to mixed states could be conceivable 
(for more recent elaborations, see [16]). 

Analysis based on the study of the dynamics of strings also support, from a different 
point of view, the idea of loss of quantum coherence. [17,18] In this respect, one can show 
that, in a rather model-independent way, a non-unitary, dissipative, completely positive 
subdynamics is the direct result of the weak interaction with an environment constituted 
by a gas of DO-branes, effectively described by a heat-bath of quanta obeying infinite 
statistics. [18] 

These results suggest that it should be possible to observe dissipative, non-standard 
effects in many quantum systems. However, crude dimensional estimates show that these 
effects are in general tiny, and therefore unobservable in practice. 

Nevertheless, a detailed analysis of the dynamics and decay of the neutral kaon sys- 
tem based on quantum dynamical semigroups shows that the non-standard contributions 
in the extended time-evolution equations could lead to testable effects. [19, 20, 21-23] These 
non-standard terms can be parametrized by six real phenomenological constants, whose 
presence modify the expressions of usual K-K observables, like decay rates and asymme- 
tries. [21, 22] Although the present experimental data are not accurate enough to detect 
these modifications, the next generation of kaon experiments should be able to put strin- 
gent limits on the six non-standard parameters. In this respect, particularly promising are 
the experiments at (/)-factories,[23] where systems of correlated neutral kaons are copiously 
produced. One should note that the description of the kaon dynamics in terms of com- 
pletely positive maps is in this case essential in order to obtain a consistent extension of 
the standard quantum mechanical time-evolution (for a complete discussion, see [24]). 



Another system in which non-standard quantum evolutions based on dynamical semi- 
groups can be studied is a neutron interferometer. [25-28] The capability of producing very 
slow neutron beams at reactors, together with the technological ability of producing and 
cutting with high precision macroscopic silicon crystals have made possible direct, very 
accurate tests of various physical phenomena. [25-30] 

In a typical experimental setup, a neutron beam is split into two components which 
travel along different paths and are subsequently brought together to interfere. The two 
components pass through a tiny slab of material before interfering; this produce a relative 
"phase shift" between the two splitted beams. By varying the relative orientation of the 
slab across the two beams, one obtains an interference figure. This figure changes under 
the action of various external phenomena, produced e.g. by earth gravity and rotation, 
or by an external magnetic field (earth magnetic field is usually properly screened); the 
corresponding theoretically calculated phase shifts induced by these phenomena have all 
been experimentally checked with high precision by analyzing the modified interference 
patterns. [28-30] 

In the following, we shall analyze in detail the dynamics of the neutron beams in 
such interferometric devices under the hypothesis that the corresponding time-evolution 
be described by a quantum dynamical semigroup. We shall see that, at least in prin- 
ciple, neutron interferometry experiments could provide very accurate estimates of the 
non-standard, dissipative terms appearing in the corresponding evolution equations. A 
preliminary analysis based on recent published data from one of those experiments will 
also be presented. 



States of a quantum system evolving in time can be suitably described by a density 
matrix p; this is a positive, hermitian operator, i.e. with positive eigenvalues, and constant 
trace. We shall analyze the evolution of neutrons in an abstract interferometer, where the 
original monoenergetic beam is splitted in two components that then interfere at the end, 
giving rise to intensity fringe patterns in two possible exit beams. We can model this 
generic physical setup by means of a two-dimensional Hilbert space, taking as basis states 
those corresponding to the two components of the split beam inside the interferometer. 

With respect to this basis, the density matrix p describing the state of our physical 
system can be written as: 

-(;:;:)■ <^' 

where ^4 = ^3, and * signifies complex conjugation. 

As explained in the introductory remarks, our analysis is based on the assumption 
that the evolution in time of the neutrons inside the interferometer is given by a quantum 
dynamical semigroup, i.e. by a completely positive, one parameter (=time) family of linear 
maps: p(0) 1— > p{t). These maps are generated by equations of the following form: 

^ = -tHpit) + ip{t)H + L[pit)]. (2) 

The first two terms in the r.h.s. of this equation are the standard quantum mechanical 
ones. They contain the effective (time-independent) hamiltonian H, that can be taken to be 



hermitian, since the fact that the neutrons are unstable can be neglected in interferometry 
experiments. 

The third piece L[p\ is a linear map, whose form is fully determined by the requirement 
of complete positivity and trace conservation: [1-3] 

L[p] = -\Y.{A^^P + pA^^) +11^opA- (3) 

J j 

The operators Aj must be such that ^ ^j^j is a well-defined 2x2 matrix; further, to 
assure entropy increase, the Aj can be taken to be hermitian. In absence of L[p\, pure 
states {i.e. states of the form |'i/')(V'|) would be transformed into pure states. Instead, the 
additional piece in (3) produces dissipation and possible transitions from pure to mixed 
states. 

As already mentioned, equations of the form (2), (3) have been used to describe various 
phenomena connected with open quantum systems; in particular, they have been applied 
to analyze the propagation and decay of the neutral kaon system. [21-24] Although the 
basic general idea behind these treatments is that quantum phenomena at Planck's length 
produce loss of phase-coherence, it should be stressed that the form (2), (3) of the evolution 
equation is independent from the microscopic mechanism responsible for the dissipative 
effects. Indeed, equations (2) and (3) are the result of very basic physical assumptions, 
like probability conservation, entropy increase, complete positivity, and therefore should 
be regarded as phenomenological in nature. 

Among the just mentioned physical requirements, complete positivity is perhaps the 
less intuitive. Indeed, it has not been enforced in previous analysis, in favor of the more 
obvious simple positivity. [10, 19, 20] Simple positivity is in fact enough to guarantee that 
the eigenvalues of the density matrix p{t) describing our system remain positive at any 
time; this is an unavoidable requirement in view of the probabilistic interpretation of p. 

Complete positivity is a stronger property, in the sense that it assures the positivity 
of the density matrix describing the states of a larger system, obtained by coupling in 
a trivial way the system under study with another arbitrary finite-dimensional one. At 
first, the requirement of complete positivity of (2) seems a mere technical complication. 
Nevertheless, it turns out to be essential in properly treating correlated systems, like the 
two neutral kaons coming from the decay of a (/)-meson; it assures the absence of unphys- 
ical effects, like the appearance of negative probabilities, that could occur for just simply 
positive dynamics. [23, 24] One should also add that standard unitary quantum mechanical 
time evolutions satisfies this property in a rather trivial way. For these reasons, in analyz- 
ing possible non-standard, dissipative effects even in simpler, non correlated systems, the 
phenomenological equations (2) and (3) should be used. 

In the particular case of neutron interferometers, as for the K-K system, a more ex- 
plicit description of (2) and (3) can be given. In the chosen basis, the effective hamiltonian 
can be written as: 

Indeed, the neutron beams inside the interferometer can be assimilated to a two- level 
system, with E the incident (kinetic) neutron energy. The splitting in energy 2uj among 



the two internal beams can be induced by various physical effects. In the following, we 
shall consider the case of a thin slab of material {e.g. alluminium) inserted transver sally to 
the two split beams. A slight rotation of this slab produces different effective interactions 
of the neutrons with the slab material in the two internal paths, yielding a non- vanishing 
uj. Using an eikonal approximation, perfectly suitable for describing slow neutrons, one 
can theoretically compute this energy splitting in terms of the neutron-nuclear scattering 
parameters for the slab material. [25-27] Typically, one finds that u is of the order of 
10"^ eV. 

The explicit form of the term L[p] in (3) can be most simply given by expanding the 
2x2 matrix p in terms of Pauli matrices ai and the identity ao: p = p^ o"^, /U = 0, 1, 2, 3. 
In this way, the map L[p\ can be represented by a symmetric 4x4 matrix [-^^jy] , acting 
on the column vector with components (po, Pi, P2, Ps)- It can be parametrized by the six 
real constants a, 6, c, a, /3, and 7: [21] 



[L,A 



/O 0\ 

a 6 c 

Q h a 13 

\0 c (3 -fj 



(5) 



(6) 



with a, a and 7 non-negative. These parameters are not all independent; the condition of 
complete positivity of the time-evolution p -^ p{t) imposes the following inequalities: 

2R = a + -f-a>0, RS > b'^ , 

2S = a + -f-a>0, RT > c^ 

2T = a + a-7>0, ST > {3 

RST >2bc/3 + RI3'^ + Sc^ + Tfe^ . 

As already observed, the dissipative correction (5) to the evolution equation (2) should 
be regarded as phenomenological; it is therefore difficult to give an apriori estimate of the 
magnitude of the non-standard parameters in (5). However, following the idea that the 
term L[p\ originates from quantum effects at Planck's scale, one expects the values of a, 
6, c, a, [3 and 7 to be very small, at most of the order m^/mp ~ 10~^^ GeV, where rrin is 
the neutron mass, while mp is the Planck scale. The dissipative contribution L[p\ in (2) 
is therefore at least three order of magnitude smaller than the one given by the standard 
hamiltonian terms: this allows an approximate analysis of the evolution equation (2), in 
which the term (5) can be treated as a small perturbation. 

For the considerations that follows, it will be sufficient to stop at the ffist order in 
the perturbative expansion. For generic initial conditions, the time dependence of the four 
components of the corresponding solution p{t) of (2) is explicitly given by: 

C C* 

piit) = (1 - 7^) Pi + ltp2 e"''^* sin(c^t) pa e*"'* sin(c^t) p^ , (7a) 

u u 

C C* 

P2{t) = 7t pi + (1 - 7t) P2 + - e-"^^ sm{ujt) ps + — e'^^ sin(c^t) p^ , (76) 

00 UJ 



P3(t) = -— e-'^' sinM) (pi - P2) + (1 - At) e-^^'^'ps + ^ sm{2ut) p^ , (7c) 

P4{t) = -- e'^' sm{ut) (pi -P2) + ^ sm{2ut) P3 + {I - At) e"'^' p^ , {Id) 

where the following convenient combinations of the non-standard parameters have been 
introduced: 

A = a + a , B = a-a + 2ib , C = c + ip . (8) 

Any physical property of the diffracted neutron beams exiting the interferometer can 
be extracted from the solution (7) for the density matrix p{t) by computing its trace 
with suitable hermitian operators. In particular, the observation of the neutron intensity 
pattern just outside the interferometer corresponds to the computation of the mean value 
of the following projector operators: [10] 

that refer to the two possible exit beams in which a neutron can be found, having traveled 
the whole interferometer; [27] the parameter ^ is a phase which depends on the specific 
experimental setup. Then, the intensity I± of the interference figure in the two exit beams 
is given by: 

I±it) = {0±) = Tr[0± Pit)] . (10) 

In the basis we are using, the initial conditions for a neutron entering the interferom- 
eter is given by one of the following two density matrices: 

They correspond to the two possible choices of orientation of the incident neutron beam 
with respect to the interferometer, and give rise to the same final results; in the following, 
we shall work with p^^\ Inserting this initial condition in the time evolution given by (7), 
from (10) one obtains the following two interference patterns: 



I±it) = l{l± 



e"^* cos {e + 2ujt) + — sin(2c^t) cos(^ - 63] 



(12) 



where \B\ and 9b are modulus and phase of B in (8); this formula holds for times such 
that At -C 1. Since a neutron having traveled fully inside the interferometer can only be 
detected in one of the two exit beams, particle conservation requires: /+(t) + I-{t) = 1, 
as it is evident from (12). 

The interference figures described by (12) are those produced by a perfectly monoener- 
getic neutron incident beam and an ideal interferometer. In practice, the neutron momenta 
in the incident beam have a finite distribution of magnitude and directions; furthermore, 
there are always slight imperfections in the construction of the interferometer, not to men- 
tion residual strains in the crystal itself. These effects can only partially be controlled, 
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and produce significant attenuation in the intensity of tfie interference figures. Detailed 
calculations based on neutron optics allow precise estimates of the modifications that are 
needed in the spectra (12) in order to take into account those effects. [25-27] In keeping with 
our phenomenological point of view, we will not use directly those estimates, but rather 
modify the expressions (12) by introducing suitable unknown parameters. By denoting 
with A^-i- the actual neutron countings at the two exit beams, one generalizes the spectra 
in (12) as: 



N±{t) 



N^f^ {1±C± 



-At 



I Rl 
cos [9 + 2ujt) + - — - sm{2ujt) cos{9 - 9b] 

2uj 



(13) 



The parameters C±, the so-called fringe contrast, take into account the previously men- 
tioned intensity attenuation, while A^^ are just normalization constants. Clearly, the 
accuracy of the determination of the parameters A and B from the measured data will 
increase as the fringe contrast gets closer to one. In actual experiments, one finds that the 
best values for C± are usually around 0.6. Further, note that now particle conservation 
requires: 



r(0) 



(0), 



N)^' C+ = N)_> C 



(14) 



In order to be able to compare the phenomenological predictions (13) with actual 
experimental data, further elaborations are required. Although the intensity spectra in 
(13) are time-dependent, in an interference experiment, being the paths followed by the 
neutrons fixed, the evolution time t is also fixed; it can only be modified by changing 
the wavelength {i.e. the energy) of the primary neutron beam. It follows that the phase 
(f = 2ujt that gives the interference figure can be varied only by changing the energy split 
u) between the two paths inside the interferometer, i.e. by changing the orientation of the 
material slab with respect to the neutron beams. 

Also, since t is fixed, it is not possible a priori to extract the fringe contrast param- 
eters C± from the (/j-dependence of the intensities A^-i-. In other words, in comparing the 
behaviour in (13) with that given by the experiment, one needs to use the following form 
for the two intensity patterns: 



iv±(vp) = ivf ai 



P± cos {9 + ip) + Q-_ 



simp 



(15) 



where 



P+ =C 



±e 



-At 



Q± = C± \B\t cos {9 - 9b) 



r(0) 



(16) 



A fit with the experimental data will give estimates for the parameters A^^ , P±, Q± and 
9; at least in principle, this is sufficient to determine the non-standard constants A and B. 
We have performed a preliminary x^ fit of the formulas in (15) with recent exper- 
imental data, published in [30]. In that experiment, a so-called skew-symmetric silicon 
interferometer and polarized neutrons were used to study the "geometrical phase" of the 
neutron wavefunction; however, in order to check the apparatus, also standard interfero- 
metric spectra had been taken. In this way, the experimental setup described in [30] is an 
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= 942 ± 6 


P+ 


= 0.17 ±0.01 


Q+ 


= 0.02 ±0.02 


^(0) 


= 366 ± 4 


P_ 


= 0.46 ±0.02 


Q- 


= 0.06 ±0.02 



actual realization of the abstract interferometer so far discussed in deriving (15). Although 
the results of other recent interferometry experiments are also available, [28] in lacking for 
an ad-hoc experiment, we find the presentation of the data in [30] the more suitable for 
the analysis of the consequences of the evolution equation (2) and (5). The results of our 
fit can be summarized as follows: 

^ = 0.09 ±0.05 (17a) 
e = 0.03 ± 0.03 . (176) 

Extracting estimates for the parameters A and B from these results requires the 
determination of the fringe contrast C±, characterizing the "optical" properties of the 
neutron interferometer. These constants are in fact function of the imaginary part of the 
refraction index of the material from which the interferometer is built. [27] Therefore, the 
best way to obtain estimates for C± is to measure the interference spectra with two neutron 
beams of wavelength A and A/2, [29] keeping unchanged the rest of the experimental setup. 
By comparing the corresponding fit estimates for the coefficients P± and Q±, with the 
help of neutron optics theory, one is then able to extract the values of C± at wavelength A. 

In lacking of a two- wavelength experiment, in the following we shall estimate C± using 
directly the data. In the standard quantum mechanical case, i.e. for A = B = 0, one 
can easily obtain the coefficients C± from the maximum A^(™^^) and the minimum A^'™'") 
neutron counts of the experimental interference figures. Indeed, from (15) with A = B = 0, 
one obtains: 

^(max) _ jY(min) 

C± = -^ — ^ f-T • (18) 

Although this relation is only approximately valid for nonvanishing A and B, in practice 
one can still use (18) with confidence, since the systematic error that one thus makes in 
the evaluation of the parameters A and B can be estimated at the end to be much smaller 
than the pure experimental uncertainty. 

Using the experimental data and (18), one obtains: C+ = 0.19 ± 0.02 and C_ = 
0.54 ± 0.03. As an independent test of the correctness of this evaluation, one can check 
that, within the errors, the relation (14) is perfectly satisfied. Note that the value of C+ 
is significally smaller than one, while that of C_ is close to the best figures that can be 
attained in practice. [26-28] This difference in the fringe contrast of the two data samples 
will result in a significally less accurate determination of the non-standard parameters A 
and B from the the results in (17a) with respects to those in (176). 

The fiight time t of the neutrons inside the interferometer can be very accurately 
determined by using the geometric specifications of the silicon interferometer used in [30] , 
so that 1/t = 5.83 x 10~^^ GeV. Further, on general grounds one expects 7le{B) and 
Im{B) to be of the same order of magnitude (for the case of the neutral kaon system, see 
[22]); then, as a working assumption, we shall neglect the small phase 9 with respect to 9b 
in (16), so that only the real part of B can be extracted from the estimates of Q±. Putting 
everything together, one finally obtains from (17a): A = (0.71 ± 0.73) x 10~^^ GeV and 
7^e(S) = (0.76 ±0.49) x lO'^i GeV, while from (176): A = (0.84 ± 0.41) x lO'^i GeV 



and TZe{B) = (0.65 ± 0.24) x 10~^^ GeV. These two estimates are compatible, but, as 
expected, the second one is much more accurate. 

Alternatively, recalling the definitions (8), one can express the previous results as an 
estimate for the parameters a and a of (5); using the best values for A and lZe{B), one 
finds: a = (0.10 ± 0.24) x lO'^^ GeV, a = (0.74 ± 0.24) x lO'^^ GeV. Although these 
values should be taken as indicative, they seem to suggest a possible nonvanishing value 
for a, while a is compatible with zero at the present level of accuracy. 

If the non-standard parameter a actually vanishes, the expression (5) of the extra term 
L[p] in the evolution equation (2) greatly simplifies. Indeed, for a = 0, the inequalities (6) 
readily imply: 7 = a, 6 = c = /3 = 0. In this case, the evolution equation (2) gives the 
most simple extension of ordinary quantum mechanics, compatible with the condition of 
complete positivity. 

In this simplified situation, the combinations A and B in (8) become both equal to a, 
so that the relations (16) are modified as: 

P± = C± e""* ~ C± (1 - at) , Q±=C± at cos 9 . (19) 

By eliminating at from these two formulas, one is now able to determine the fringe contrast 
factors C± from the fitted parameters P± and Q± without further assumptions: C+ = 0.20± 
0.02, C- = 0.52 dz 0.03. Note that these values are equal within errors to those determined 
before. Then, from either of the relations in (19), one obtains two determinations of the 
parameter a, which combined finally give: 

a = (0.71 ±0.21) X 10-21 GeV . (20) 

Although this estimate for a points toward a nonvanishing value, roughly of the right 
order of magnitude for a quantum gravity or "stringy" origin, it should not be regarded 
as an evidence for non-standard, dissipative effects in the dynamics describing neutron 
interferometry. Rather, it should be considered as a rough evaluation of the sensitivity 
that present neutron interferometry experiments can reach in testing quantum dynamical 
time evolutions of the form given in (2) and (5). 

In closing, we would like to make a few comments on the existing literature on the 
subject. In our study of the effects of the environment on the propagation of the neutrons 
inside the interferometer, we have assumed that the refractive phenomena on the various 
silicon blades of the device be described by standard neutron optics. This theory is the 
result of a quantum mechanical analysis of the scattering of the neutron beams by the 
nuclei in the silicon crystal. For slow neutrons, the effects of these scattering phenomena 
can be effectively described by a model that resemble very closely "geometric optics" of 
light propagation theory. [25-27] 

In principle, non-standard dissipative effects can also be present in the scattering 
neutron-nucleus; [31-33] these phenomena can be described again by phenomenological 
equations of the form (2), (3), and could modify the predictions of standard neutron optics. 
A precise estimate of these changements requires detailed computations that certainly go 
beyond the scope of the present investigation. In any case, it should be stressed that these 
possible dissipative effects in the interaction neutron-nucleus would mostly affect the es- 
timate of the intrinsic interferometer parameters, like the fringe contrast or the phase 9, 
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as functions of the wavelength of the incident neutrons and the nuclear properties of the 
refractive material. In our study, these parameters have been obtained directly from the 
experimental data, so that we expect little changements in our analysis from these extra 
effects. Nevertheless, the entire topic certainly deserves further attention and we hope to 
come back to these problems in the future. 

A study of possible phenomena violating quantum mechanics in neutron interferome- 
try has been originally presented in [10]. There, an equation of the form (2) has also been 
used to describe these effects, but without imposing the condition of complete positivity. 
Fitting an approximated formula for the exit beams interference figures with the experi- 
mental data available at that time, limits on some of the parameters violating quantum 
mechanics were given. These limits have been further strengthen by later analysis, [34] 
exploiting wavefront splitting interference experiments (the analogs of Young's two slit 
experiment). 

These estimates turn out to be rather crude: they are based on a rough evaluation 
of the flight-time of the neutrons inside the interferometer, rather than a detailed analysis 
of the interference patterns. Furthermore, as already mentioned, lacking of imposing the 
condition of complete positivity on the evolution equation could lead to serious inconsis- 
tencies. [24] We stress that to avoid these problems, one needs to adopt phenomenological 
descriptions based on equations (2) and (3). 

Finally, the neutron interference experiments realized so far allows determining at 
best the values of only two of the six non-standard parameters in (5). The remaining ones 
could be estimated, at least in principle, by studying the behaviour of other observables 
O, different from those appearing in (9). This would allow a direct test of the inequalities 
(6) and therefore of the hypothesis of complete positivity. 

In practice, however, the analysis of these new observables would correspond to the 
realization of completely different experimental setups. In this respect, a detailed study 
of possible non-standard, dissipative effects in neutron interferometry appears to be an 
exiting challenge not only theoretically, but experimentally as well. 
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